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Abstract 

Little  or  no  test  data  is  available  to  verify 
theory  predicting  lift,  drag  and  pitching  moment 
of  a  hydrofoil  in  the  presence  of  pitching  motion 
coitmonly  encountered  in  hydrofoil  craft.  The 
present  paper  presents  the  results  of  a  wind  tun¬ 
nel  investigation  of  the  effects  of  pitching  motion 
on  two  hydrofoil  configurations.  The  results  are 
compared  to  prediction  techniques  given  in  the 
U.S.A.F.  OATCOM  (Ref.  1).  Finally,  a  verified 
force  and  moment  prediction  technique  is  presented. 


1 .  Introduction 

A  great  deal  of  emphasis  is  placed  on  the 
Ditching  motions  of  conventional  displacement 
ships,  i.e.,  on  methods  of  prediction,  and  methods 
of  reducing  the  severity  of  these  motions.  These 
same  motions  occur  in  no  less  a  degree  on  high 
performance  hydrofoil  craft,  in  this  case  however 
pitching  motion  is  encountered  by  the  hydrofoils, 
rather  than  the  hull  proper.  Unfortunately,  in 
dynamic  simulations  of  these  craft,  the  effects  of 
pitching  motion  are  notably  absent.  (References 
2  and  3)  This  is  due  to  the  lack  of  an  experi¬ 
mentally  verified  theory  on  the  hydrodynamic 
properties  of  hydrofoils  in  pitching  motion.  Such 
a  theory  would  have  wide  ranges  of  application  in 
stability,  control  and  seakeeping  studies  of  hy¬ 
drofoil  craft.  Further  areas  of  application  im¬ 
mediately  come  to  light,  including  studies  on 
ship  rudders  in  yawing  motion,  submarine  control 
surfaces  in  both  pitching  and  yawing  motions  and 
also  in  S.W..A.T.H.  ship  stabilizer  and  strut 
surfaces  in  both  pitch  and  yaw. 

A  survey  of  the  literature  indicates  rela¬ 
tively  little  work  in  this  area.  Ribner  (Ref.  4) 
gives  approximate  methods  of  obtaining  both  static 
and  dynamic  derivatives,  but  these  are  limited  to 
triangular  planforms  with  low  aspect  ratios.  Toll 
and  Queijo  (Ref.  5)  presented  the  approximate 
method  of  considering  the  hydrofoil  lift  and  mo¬ 
ments  to  be  the  sum  of  three  components, 

1.  an  increment  of  lift  and  pitch  due  to 
the  aerodynamic  center  being  offset 
from  the  center  of  gravity,  about  which 
pitching  motion  is  defined, 

2.  an  increment  of  lift  and  pitch  due  to 
the  curvature,  approximated  by  assuming 
the  curvature  to  cause  the  foil  to  be¬ 
have  as  if  it  had  a  circular  camber,  and 
assuming  the  additional  lift  to  act 
through  the  midcord, 

3.  an  Increment  of  lift  and  pi  ten  due  to  the 
additional  lift  of  tne  above  increment, 
acting  throuon  the  quarter  chord,  but  in 
the  opposite  direction. 
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This  method  is  presented  for  use  in  the  DATCOM 
(Ref.  1).  Calculated  dynamic  derivatives  using 
this  method  are  compared  to  the  data  from  the  cur¬ 
rent  investigation. 


The  present  paper  gives  the  results  of  curved 
flow  wind  tunnel  tests  on  one  low  aspect  ratio 
and  one  moderate  aspect  ratio  hydrofoil.  These 
data  are  then  used  in  a  proposed  semi -empirical 
model  of  hydrofoil  properties  under  the  Influence 
of  pitch  rates. 


The  investigation  was  conducted  in  the  V.P.I. 

&  S.U.  subsonic  stability  wind  tunnel  using  the 
6'  x  6'  curved  flow  test  section.  In  this  sec¬ 
tion,  the  vertical  tunnel  walls  are  physically 
curved  to  provide  flow  angularity  while  wire 
screens  creating  a  dynamic  pressure  gradient  across 
the  test  section  provide  the  required  velocity  dis¬ 
tribution  to  complete  the  simulation  of  curved 
motion.  Although  the  test  section  walls  may  be 
formed  to  a  wide  range  of  curvatures  the  number  of 
wire  screens  available  limit  the  achievable  turn 
rates  to  four.  Defining  a  nondimensional  pitch 
••ate  as  q'  *_qc/2 V,  where  q  is  the  dimensional 
pitch  rate,  c  is  the  hydrofoil  mean  aerodynamic 
chord  and  V  is  the  freestream  velocity,  pitch  rates 
of  0.0387,  0.0572,  0.0787,  and  0.1115  (based  on  a 
chord  of  4  feet)  can  be  attained.  A  more  detailed 
description  of  the  facility  together  with  its 
theory  of  operation  is  given  in  Reference  6. 


Although  the  test  section  is  constructed  to 
give  yaw  rates,  pitching  may  be  simulated  by  ro¬ 
tating  the  hydrofoil  90  degrees  to  an  upright 
position.  In  the  present  investigation,  two  hydro¬ 
foils  were  tested  in  two  curvatures  numbers  two 
and  four,  and  straight  flow.  Table  1  gives  the 
various  physical  particulars  of  each  foil  to¬ 
gether  with  tha  pitch  rates  corresponding  to  the 
tunnel  curvatures  at  which  the  foils  were  tested. 
The  foils  were  mounted  vertically  on  a  strut  in 
the  tunnel  on  top  of  a  6'  x  4'  ground  board  with 
approximately  .25  Inches  clearance  above  the  board. 
The  ground  board  serves  as  an  endplate  to  ef¬ 
fectively  double  the  aspect  ratios  of  the  models. 

At  each  curvature,  the  models  were  rotated  up  to 
angles  of  attack  of  plus  or  minus  10  degrees  in 
increments  of  two  degrees.  Mounted  on  the  strut 
Inside  each  foil  was  a  six  comoonent  strain  gauge. 
At  each  angle  of  attack  forces  and  moments  were 
read  and  reduced  using  a  Hewlett  Packard  Data 
Acquisition  System. 


3.  Results 

In  the  present  work,  the  particular  interest 
is  on  the  behavior  of  lift,  drag  and  pitching  mo¬ 
ment  in  curved  flow.  These  quantities  were  reduced 


from  strain  gauge  measurements  by  Including  calcu¬ 
lation  of  blockage  effects  due  to  the  presence  of 
the  large  hydrofoils  (compared  to  tunnel  cross 
sectional  area)  and  the  calculation  of  bouyancy 
terms  due  to  the  presence  of  the  radial  pressure 
gradient  caused  by  the  tunnel  curvature.  These 
forces  and  moments  were  nondimensional ized  by  the 
following  conventions i 


L/(l/2)pV2S 

(1) 

D/(1/2)pV2S 

(2) 

m/(l/2)pV2S 

(3) 

where  l,  D,  and  M  are  the  lift,  drag  and  pitching 
moment,  respectively;  S  is  the  planform  area,  and 
Cl,  Cq,  and  Cm  are  the  nondimensional  lift,  drag 
and  pitching  moment.  Table  1  lists  the  location 
of  the  reference  center  about  which  the  pitch¬ 
ing  moment  was  measured. 


TABLE  1. 


Foil  1 

Foil  2 

Section 

MAC* ** ***  0020 

NACA  0020 

E9HI 

2.17ft. (,83m) 

1.4ft. ( .34m) 

c* 

1 .B6ft.(  .57n) 

2.37ft. ( .h5m) 

Taper  ratio 

0.5 

1.0 

AR 

2.9 

1.45 

tmm 

l. 125ft. (.35n) 

C.75ft.( ,23n) 

Sweep  of  :<c 

20.34  (deg. ) 

0.0  (dee.) 

*2 

0.0266 

0.0339 

u 

0.052 

0.066 

*  Mean  Aerodynamic  Chord 

**  01  stance  from  foil  Vertex  to  moment  reference 
center 

***  Nondimensional  pitch  rate  at  curvature  2  and  « 


Figure  l.  Lift  Coefficient  vs.  Annie  of  AtUck  For 
Foil  1. 


DATCOM  (Reference  6).  These  comparisons  are  given 
in  Table  2.  Examination  of  these  shows  good  agree¬ 
ment  of  V.P.I.  calculated  CLa  and  Cmq  with  the 
DATCOM  predictions.  The  agreement  of  V.P.I.  cal¬ 
culated  Cm  and  Ctq  with  DATCOM  methods  is  seen 
to  be  relatively  poor.  This  is  probably  due  to  a 
difficulty  in  calculating  the  location  of  the  aero¬ 
dynamic  center. 


Figures  1  and  2  show  plots  of  Cl  versus  angle 
of  attack  for  foils  1  and  2  for  all  three  cur¬ 
vatures  tested.  The  nonzero  intercept  of  the 
zero  curvature  line  is  due  to  the  presence  of 
some  angle  of  attack  due  to  misalignment  of  the 
foils  in  the  tunnel.  In  subsequent  data  analysis 
these  misalignments  were  treated  as  tares. 

Figures  3  and  4  are  plots  of  Cm  ve>-«ns  angle 
of  attack  for  foils  1  and  2  showing  al1  three 
curvatures.  Foil  misalignment  is  ind-'  jted  by  a 
repeated  upward  shift  of  the  zero  cur*,  ‘ure  data. 

From  this  data,  stability  derivatives  were 
calculated  and  compared  to  those  predicted  by 
the  method  of  Toll  and  Quel jo  as  given  in  the 


TAM.[  Z. 


V.P.I.  1  s.u. 

DATCOT' 

foil  1 

Foil  2 

Foil  1 

Foil  2 

C,  ,(/r«4) 

3.34 

I.SO 

3.41 

1.S7 

C^.(/rid) 

0.54 

o.s« 

1.23 

0.32 

C,.(/r»4) 

0.0 

0.0 

-1.436 

-0.245 

c_(/r«a) 

-0.75 

-0.S4 

-0.7Z 

•0.62 

Figure  2. 


lift  Coefficient  vs.  Annie  of  Attack 
for  foil  2 


Figure  1.  fitch  Moment  Coefficient  es.  Angle  of  Attack 
for  foil  1. 


4.  Modeling 

4.1  Introduction 

Currently,  in  linear  dynamic  simulations,  the 
vehicle  forces  and  moments  are  expanded  using  a 
Taylor  series  and  taking  the  motions,  accelera¬ 
tions  and  displacements  to  be  perturbation 
variables  (Ref.  7).  Terms  of  order  greater  than 
one  are  generally  neglected  in  the  Taylor  series. 

In  an  expansion  of  the  lift  coefficient,  the  as¬ 
sumed  form  may  be,  *  Ci_(a,q',r', . . . ) ,  and  the  cor¬ 
responding  Taylor  expansion,  Cl  *  Cio  +  3Cl/ 

3a*i  +  3CL/aq,*ql+  ...  etc.,  where  is  the 

3a  o 

stability  derivative.  The  first  two  terms  of 
this  expansion  are  generally  well  understood  and 
have  a  very  clear  physical  Interpretation.  How¬ 
ever,  the  third  term  is  somewhat  less  well  under¬ 
stood.  Essentially  the  effects  of  pitch  rate  are 
grouped  together  in  one  term  (Ref.  8)  and  ft  is 
difficult  to  gain  insight  into  the  origins  of  the 
basic  hydrodynamics.  The  intent  of  the  proposed 
models  is  to  allow  description  of  the  effects  of 
pitching  motion  on  hydrofoils  in  a  manner  that 
arises  more  naturally  from  the  data  and  also  al¬ 
lows  these  effects  to  be  expressed  in  terms  of 
explicit  parameters. 


4.2  Lift  Coefficient  Development 

To  model  the  lift  coefficient  to  include 
pitching  effects,  the  first  two  coefficients  of 
the  Taylor  expansion  in  angle  of  attack  were  as¬ 
sumed  to  be  arbitrary  functions  of  q'.  This  al¬ 
lows  the  camber  inducing  effect  of  pitching  as 
suggested  by  Toll  and  Queijo  (Ref.  5),  as  well  as 
the  behavior  of  the  conventional  lift  curve  slope, 
to  be  seen.  The  model  describing  the  lift  co¬ 
efficient  linear  in  a,  is: 

CL  -  CL  +  aCL  (4) 


where,  in  general,  the  Intercepts  Clq  and  the 
slopes  CLa  are  functions  of  the  steady  pitching 
rate  q'.  From  figures  1  and  2  it  is  apparent  that 
CLp  I*  very  nearly  independent  of  q'  and  Cl„  is 
likewise  constant  with  changine  q'.  With  both 
coefficients  of  equation  4  constant,  the  derivative 
with  respect  to  q'  becomes: 

C.  •  ■  3C./30'  *  0  (5) 

Lq  L  a*0 


of  pitch  rate.  Therefore  the  model  describing 
drag  behavior  In  pitching  motion  becomes: 


Figure  4.  Pitch  Moment  Coefficient  Angle  of 
Attack  for  foil  2 


4.3  Drag  Coefficient  Development 

In  modeling  the  drag  coefficient,  the  drag 
polar  equation  was  used,  allowing  the  drag  to  be 
expressed  as  a  function  of  zero  lift  drag  and 
lift  due  to  Induced  drag.  The  following  equation 
was  then  used: 

CD  *  C0o  +  7k  <Cl  *  Cl  <6> 

O  Oj 

where.  In  general,  Cq0.  Cl0  and  k  are  functions 
of  the  pitch  rate  q'and  AR  Is  the  aspect  ratio. 

Figures  S  and  6  are  plots  of  drag  polars  of 
each  foil.  These  plots  allow  Isolation  of  the 
parameters,  C|n0  and  CL02  that  vary  with  q'.  Fur¬ 
ther  manipulation  of  the  data  to  Identify  a 
variation  of  the  parameter  k  with  q'  led  to 
plotlng  the  drag  polar  In  the  form;  (Co  -  CDo ) 
against  (C|.  -  Cl?).  This  form  allows  k  to  appear 
as  the  slope  of  the  now  linear  curve.  As  seen 
In  figures  7  and  8,  k  turns  out  to  be  Independent 


co*co<<,)  +  ;*&<cL-ci<<’,»2  <7> 

0  o 

and  the  corresponding  derivative  with  respect  to 
q'  becomes: 

3C0/3q'  *  ®C0Q/3q '  *  *Cl  *  CL0}  *CL0/#q' 

(8) 

Figure  9  shows  Cl0  and  Cp0  as  functions  of  q'  for 
each  foil.  Using  this  plot,  the  derivatives 
3C0o/39'  and  3Ci0/3q'  may  be  calculated  and  em¬ 
ployed  In  equation  7  to  find  sCo/sq'.  It  Is  seen 
that  the  dependence  of  Cl0.  Cd0  and  k  are  strongly 
Influenced  by  physical  parameters  of  the  foils 
such  as  aspect  ratio,  taper  ratio,  thickness,  etc. 
Unfortunately  Insufficient  data  was  collected  to 
completely  characterize  this  Influence.  The  col¬ 
lection  of  data  in  which  these  parameters  are 
varied  systematically  would  generalize  the  present 
approach  and  allow  an  extended  semi-empirical 
method  to  be  developed. 


4.4  Pitch  Moment  Coefficient  Development 

The  pitching  moment  was  modeled  using  an 
equation  of  a  form  analogous  to  that  used  in  the 
lift  coefficient  model.  Therefore  the  following 


o  Curwtur*  0 


n«nrt  t.  llftMrlittf  Orif  Ptlir  for  foil  l. 


M<hir*  y.  Linoartifd  >*f  *ol»f  »0f  #0l1  1. 


equation,  linear  with  angle  of  attack  was  used: 


+  C_ 


(9) 


where,  C,^  and  are  functions  of  q'. 


flforo  f.  Cqq  ooO  C^#  for  f»11»  1  t. 


Figures  3  and  4  indicate  that  both  of  these 
parameters  are  dependent  on  q\  This  dependence 
Is  shown  for  the  two  foils  in  figure  10.  The 
model  describing  pitching  moment  behavior  In  pitch, 
ing  motion  becomes: 

Cm“CmB(q')+C«(,,,,‘  U0) 

o  a 

and  the  derivative  *  Is: 

VV*,Ch/Sq'MC«  'V*  (ID 

0  o  ^ 


2 


As  In  the  development  of  the  drag  coefficient, 
the  terms  in  equation  10  may  be  calculated  from 
the  plot  of  the  respective  functions  (Figure  10). 
Again  it  is  seen  that  both  and  Cm,  are 
strong  functions  also  of  the  physical  parameters 
of  the  foils. 


Flrure  10.  C  <nd  C„  »or  Foil!  1  «n<l  3. 
no 
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5.  Conclusions 

Hind  tunnel  data  of  hydrofoils  in  curved 
flow,  simulating  pitching  motion,  has  been  used 
in  developing  new  models  for  describing  the  ef¬ 
fects  of  pitch  rate  on  hydrofoils.  Expressions 
for  lift,  drag,  and  pitch  moment  coefficient  have 
been  presented  which  Isolate  the  Influence  of 
pitch  rate  on  parameters  more  commonly  encoun¬ 
tered  and  more  physically  tractable  in  hydrody¬ 
namic  work.  It  is  felt  that  continued  investi¬ 
gations  in  which  physical  characteristics  of  the 
hydrofoils  are  varied  would  produce  very  useful 
extensions  of  the  proposed  models. 
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